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SUMMARY 

The first two terms in the general spherical harmonics expansion of the 
Fokker-Planck equation are derived hy a simplified method. Some remarks are 
made on the feasibility of obtaining the general term of the expansion. Errors 
in previously published results are pointed out. 


INTRODUCTION 

The Fokker-Planck equation is extensively used in the analysis of non- 
equilibrium phenomena in plasmas. The basic assumption underlying this equa- 
tion is that collisions which result in scattering through small angles are the 
most probable and therefore have the predominant effect on the distribution 
function. This assumption is well founded for electrically charged particles 
interacting according to an inverse -square force. 

If the validity of the Fokker-Planck equation is accepted, a difficult non- 
linear integro -partial differential equation for the distribution function still 
remains. The usual method of simplifying it involves the expansion of the dis- 
tribution function in spherical harmonics in velocity space. The result of 
this procedure is a system of coupled nonlinear integro-differential equations 
whose independent variables are the magnitude of the velocity and time. 

Rosenbluth, et al. (ref. l) described the spherical harmonics expansion 
for the case of axial symmetry, but the resulting system of equations was not 
presented. In a subsequent paper (ref. 2), however, the equation for an iso- 
tropic distribution function, which is the zero-order equation of the system, 
was given. Dreicer (ref. 3) obtained the zero- and first-order equations of 
the system. Other authors (refs. 4 and 5) have also derived the zero- and 
first-order equations. These two equations usually suffice for most practical 
purposes, since the zero-order term in the spherical harmonics expansion gives 
the energy distribution of the particles, while the first-order term describes 
the manner in which the current is distributed over particles of various ener- 
gies. 


The Fokker-Planck equation for electrons may be combined with the usual 



electric and magnetic field terms of the Boltzmann equation and the appropriate 
collision terms for elastic and inelastic interactions Between electrons and 
molecules to yield a description of phenomena in a partially ionized gas 
(ref, 3). Even in a partially ionized gas, electron-electron collisions have a 
significant effect on the electron distribution function, except when the degree 
of ionization is extremely low. The inclusion of this effect entails consider- 
able analytical difficulty unless simplifications of the type described previ- 
ously can be made. Because of the importance of electron-electron collisions 
in many situations, the correct formulation of the collision term is vital. 

Errors of a nontrivial nature were discovered by the present author in the 
results of references 3 to 5. Dreicer T s results (ref. 3) are the most nearly 
correct, but much more serious errors occur in references 4 and 5. The ana- 
lysis leading to these results is not simple, but it could not be included in 
detail in the previously cited papers because of space limitations. Therefore, 
it is worthwhile to present in full the derivation of the first two equations 
in the spherical harmonics expansion of the Fokker-Planck equation and to point 
out the errors in the results of references 3 to 5. Although the derivation 
follows the general outline of references 1 and 3, parts of the derivation 
presented herein are believed to be somewhat simpler than the analyses described 
in those references. 


FOKKER-PLANCK EQUATION 

The Fokker-Planck equation may be written (refs. 1 and 3) 


'Sf\ 


((AVjJF) + | ((AViAYj)*) 


m l 


(la) 


where F(v,r,t) is the distribution function, which is defined such that Fd 3 v 
is the number of test particles in the element d^v in velocity space centered 
about the velocity v. The components of v on a fixed system of rectangular 
Cartesian coordinate axes are denoted by v^(i = 1, 2, 3). Repeated subscripts 
imply summation. (Symbols are defined in appendix A.) 


Unfortunately, the terminology applied to the various kinetic equations is 
by no means universal. Equation (la) will be referred to as the Fokker-Planck 
equation in this report, following the usage of references 1 and 3. Actually, 
equation (la) is simply the definition of (c)F/St) c , the rate of change of F 
due to collisions. In the general case where there are spatial gradients and 
external forces acting on the test particles, expression (la) is equated to the 
total derivative following the motion of the test particles 


5f 
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and equation (lb) is sometimes called the Boltzmann-Fokker-Planck equation. 

The vector (AVj_) is the rate of change in the i^ h component of velocity of 
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k 



Figure 1. - Orientations of velocities 
of test and field particles with re- 
spect to coordinate system. 


a test particle of velocity v-^ due to its simultaneous 
interaction with all other particles, which are called 
field particles. If the distribution function of the 
field particles is denoted hy F f , (Av^) is defined as 
(refs. 1 and 3) 


<avi>= Zi v- ,)dV x 


o{g,X)g Av ± d£2 


( 2 ) 


where 

scattering through the angle 


cr(g,X) is the differential cross section for 
X, g is the relative speed, that is, 


g 



(3) 


Av^ is the increment in velocity of a test particle of initial velocity v 
due to a collision with a field particle of velocity v T , and da is the ele- 
ment of scattering solid angle. The velocity vectors v, v ! , and g are 
shown in figure 1. The summation in equation (2) is carried out over all 
species of field particles. The meaning of equation (2) becomes clearer if it 


is noted that 



f 



a g da 


is the instantaneous collision rate 


of a test particle of velocity v with all of the field particles. The quan- 
tity (Av-j_) is frequently called the coefficient of dynamical friction, since 
it corresponds to the slowing down of the test particles along their direction 
of motion (ref. 6). Note also that (Av_^) has units of velocity per unit time. 


The quantity (AvqAv.) is defined similarly as 

J 


( ^> .y i F f (v , )d 3 v ! j* a(g,X)g Av^Avj 


da 


( 4 ) 


and is referred to as the dispersion coefficient, since only this term in the 
Fokker-Planck equation contributes to the spreading of a beam of test particles 
(ref. 6). 

The Fokker-Planck equation (eq. (la)) may be derived from the Boltzmann 
equation by means of an expansion of the kernel of the Boltzmann collision 
integral in powers of the velocity increments (ref. 6). The Fokker-Planck equa- 
tion is obtained when only terms through the second order in the velocity in- 
crements are retained. Thus the Fokker-Planck equation is valid when collisions 
resulting in small angle scattering are predominant. This is true for Coulomb 
interactions because of the long range nature of the inverse -square force. 
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The differential cross section for Coulomb scattering is the familiar 
Rutherford cross section 


o( g,X) 




In this equation the' subscripts t and f refer to the test and field parti- 
cles, respectively. The reduced mass m^ is 


m = 




tf m^. + m^ 


( 6 ) 


(Mks units are used throughout this report. Thus the electric charges e^_ 
and e^ are measured in coulombs (e.g. , e^ = ±Z^e where is the 

charge number and e = 1. 602xl0~ x C is the absolute value of the electron 
charge). All the results herein may be converted to cgs electrostatic units 
by simply deleting the 4 k€ 0 .) 

The calculation of (Av^) and (Av^AVj) for Coulomb collisions was carried 
out by Rosenbluth, et al. (ref. 1), who showed that these quantities could be 
expressed by means of the following simple formulas: 

< Av i> - < 2 * * * * 7 > 

i 



dv.dv . 


( 8 ) 


where the functions H^. and G^, which are now generally known as the 
Rosenbluth potentials, are defined as^ 


Treicer’s expression for a (following eg. (4) of ref. 3) erroneously has 

. o o,2 

in the numerator. This rather trivial mistake is not one of the errors 

"G I 

referred to in the INTRODUCTION. 

2 

Dreicer’s notation is used in equations (9) to (ll) to facilitate com- 

parison with his final results. His definition of and differs from 

that of Rosenbluth, et al. (ref. l) by the presence of r tf Dreicer ! s equa- 

tion for r tf (eq. (8) of ref. 3) contains the same error as noted in foot- 
note 1, in addition to incorrectly having the reduced mass instead of m-^. The 

latter error also appears in reference 1. 
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( 9 ) 


H + (v) = 


f 


m f ^F f (v')d 3 v’ 

r tf / g 


G t (v) = G f (v) = T tf j* F f (v')g <a 3 v' 
f f 

and the velocity- independent factor is 


( 10 ) 


r« - - ( S ?) 2 ^ " - ( S ?) 2 ^ 

t t 


In A, 


( 11 ) 


The quantity A^, which is usually much larger than unity, is defined as 
the dimensionless ratio 


A tf - 


( 12 ) 


•where d is the Debye radius and b Q the impact parameter for 90° scattering. 

When the charged species are in thermodynamic equilibrium at the temperature T, 
the Debye radius is given by (ref. 6, p. 2l) 



(13) 


where the summation is over all charged species. The average impact parameter 
for 90° scattering is (ref. 6, pp. 16 and 21) 


b 


o 


e t e fi 

4jre„ 



e t e fl 

4:rre 


1 

3kT 


(14) 


Q 

since for thermodynamic equilibrium, m^(g^) = 3kT, as may be easily verified 
by integrating g2 over IViaxwe Ilian distributions of field and test particles. 
Therefore, equation (12) becomes 
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By direct differentiation of equations (9) and (10), the following useful 
relations can he derived 3 (refs. 1 and 3) 


?vH t 


a 2 H t 

Sv ;i _dv i 



r tfV 7 > 


7 


2 


v 



(16) 


(17) 


-,4= n - ^ 4 ^t 

V V ^t c)vj_dvj_dvj dvj 



(18) 


where 


Hp and H, are related hy the first identity in equation (9). 

X U 


With the use of equations (7) and (8), the Fokker-Planck equation 
(eq. (la)) for Coulomb collisions may be written in terms of the Rosenbluth po- 
tentials H, and G, as follows: 



dv . dv • 


a2s t \ 

dT ± d d/ 


(19) 


Rosenbluth (ref. l) and Dreicer (ref. 3) express equation (19) in spherical 
coordinates without further reduction. For the present purposes , it is more 
convenient to first reduce equation (19) to a simpler form. Expansion of the 
derivatives of the products leads to 



r SF SH t + 1 p s4q t x dF 3 / s2& t V 1 6 2 F 52& t / 2Q) 

SvqSvq dvj_ cFvq 2 Sv^SvT^VjBv^ dvj_ \dvjdvjj 2 dv-jdvj dvqdvj ' ' 


3 The factor r._ is missing from equations (9a) and (9c) of reference 3. 

"Cl 
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Substitution of equations (16) to (18) into equation (20) yields 


3F S - m f „ , 1 8 2 f a2 °t 

( v ) 5v7 dv7 / m , + nu f 2 dv.dv. 5v Bv 

i 1 / -j t ± 1 J i j 

f f (21) 

This form of the Fokker-Planck collision term is not ne¥ ; since it was given by 
Longmire (ref. 7), for example. Equation (21) is converted to an equation in 
spherical coordinates most simply by first writing it in terms of general co- 
variant derivatives valid for any coordinate system and then specializing the 
result to spherical coordinates. This procedure is carried out in appendix B. 
The result, equation (B17), could be used as the starting point of a general 
spherical harmonic expansion of the Fokker-Planck equation. Because of the 
nonlinear form of the equation, however, the general term of the spherical 
harmonics expansion would be exceedingly cumbersome and not too useful for 
practical purposes. This report will therefore confine itself to the deriva- 
tion of only the first two terms of the spherical harmonics expansion, and for 
this limited purpose a much less cumbersome procedure will suffice. 


= 4rt 


cc 


i 


“t 

— r f 

m f tff 


EXPANSION IN SPHERICAL HARMONICS AND THE LORENTZ APPROXIMATION 


A considerable simplification of the Fokker-Planck equation is obtained if 
it is assumed that the distribution function is nearly spherically symmetric in 
velocity space. The distribution function is written as a spherically symme- 
tric distribution F°(v) (not necessarily Maxwellian) plus a small perturbation 
i(v); that is 

F(v) = F°( v) + l(v) (22) 


This is equivalent to retaining the first few terms in the general spherical 
harmonic expansion of F 

F(v) ^ F 2m (v)Y Zffi (e,cp) (23) 

Z= 0 m=-l 

In general, F^ m is a function of the position vector r in physical space and 
time, as well as v. This space and time dependence is understood in equa- 
tion (23) and all subsequent expressions. The surface harmonics Y, (0,cp) are 
defined as (ref. 8, sec. 3.4) 

W'O - J zl l ^ 0 > elKp <*•> 


where are the associated Legendre polynomials , defined as 

L 


_m, ^ 1 , nl m d l+m (cos 2 0 - 1) 

P 7 (c°s e) = 7Y77 (- sin 0 ) 7T w.+ 5T 


2 b ll 


(d cos 0) 


(24b) 
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The angles 6 and cp are the polar and azimuthal angles, respectively, mea- 
sured with respect to an arbitrary hut fixed orthogonal triad of unit vectors 

i, j, k (fig. l). Since = (-l) 11 ^^ where th e asterisk denotes the com- 

plex conjugate, one way to ensure that F is real is to assume that 
F Z,-m = ( _ l) m (F Zm ) * 

Although the expansion of F (eg., (23)) is convenient for many purposes 
because of the orthogonality properties of the spherical harmonics, it still 
must be interpreted with reference to some set of axes, and this is sometimes 
inconvenient. It is often desirable to recast equation (23) in an invariant 
form, which is independent of the particular choice of axes. This was done by 
Johnston (ref. 9), who has shown that the expansion of F (eg. (23)) is equiva- 
lent to a single sum of contracted tensors of the form 


00 



z=o 


where F. . . is a tensor of order Z. The general expression relating 

Xli 2 . . .lj 

this tensor and the coefficients F^ m of equation (23) is given in reference 9. 
This somewhat complicated expression will not he required in the present report 
because only the first two terms ( Z = 0, l) of equation (23) will he considered. 
When these terms are written out explicitly hy using the expressions for Y^ m 
obtained from equations (24), 


1 



Equation (26) may he cast into the desired form 

F(v) = f®(v) + Fl(v) • — + . . . (27) 

v 
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(28a) 


by making the definitions 4 

F 00 (v) 

V2 ImjF i;L (v)J j + F 10 (v)k | (28h) 

When there is rotational symmetry about the polar axis £, only the terms 
for m = 0 in equation (23) are retained and the expansion reduces to an ex- 
pansion in ordinary Legendre polynomials 



F(v) = F Z (v)P z (cos 9) 


where 


F Z (v) = | / 2Z 4 * 1 F Z0 (v) 


(29) 


(30) 


(see eq. (3.57) of ref. 8). In this case, the vector F x is directed along 
the polar axis, as is clear from equation (28b). For axial symmetry, the polar 
axis is the direction of the force on the test particles, for example, an 
electric field or concentration gradient (see ref. 10, pp. 404-408). In the 
presence of electric and magnetic fields, however, there is in general no 
axial symmetry in velocity space because, for arbitrary directions of 

E and B, the force (E + v X B) is not rotationally symmetric about any axis 
and the direction of F (v) is unknown a priori. 


The Lorentz approximation (see refs. 3 and 11) consists of retaining only 
the first two terms in equation (27) under the assumptions that. 



« F° 


(31) 


and that the remaining terms are of negligible importance. As Johnston (ref. 9) 
has pointed out, F^ alone is needed to calculate the average value of any 
scalar function S(v), such as the kinetic energy, while the average value of 
any vector quantity of the form Q(v)v, such as the velocity vector itself or 

-»1 

the heat flux vector, depends only upon F . Terms of higher order than those 
written out in equation (27) account for higher order tensor properties, such 
as the nondiagonal elements of the pressure tensor or the energy flux tensor 
(ref. 9). 


The physical meaning of condition (31) can be established as follows. The 
average drift velocity of the test particles is (ref. 9, eq. (10b)) 


^Equations (28a) and (28b) differ from equations (9a) and (9b) of refer- 
ence 9 by numerical factors because the spherical harmonics are defined dif- 
ferently in reference 9. 
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v 3 F 1 (v)dv 


(32) 





CO 


The average random speed (v) is (ref, 9, eq. (10a)) 


<v) 



(33) 


Therefore, condition (3l) implies that the magnitude of the drift velocity is 
much smaller than the random speed. 

It should he emphasized that in the present analysis, as well as in the 
analyses of references 1 to 3, neither the spherically symmetric part F^(v) 
of the distribution function of the test particles, nor the corresponding func- 
tion F^(v) for the field particles, is assumed to be Maxwellian. Furthermore, 

the effect of collisions among the test particles ( self -collisions ) is included 
in the analysis, inasmuch as the summations over the various species f of the 
field particles include the case f = t, where the field and test particles are 
the same. In general, the Fokker-Planck equation is nonlinear due to the pres- 
ence of the self -collision terms. 

On the other hand, some authors (refs. 12 and 13) linearize the Fokker- 
Planck equation by assuming that the distribution functions of both the test 
and field particles may be represented by a Maxwellian distribution plus a 
small perturbation. In the notation of this report, such an assumption cor- 
responds to taking F°(v) and F^(v) to be Maxwellian distributions, for which 
the rate of change of fO(v) due to collisions vanishes. Linear ordinary 
integro-differential equations for the perturbation terms, such as F^-(v) (in 
the notation of this report), are obtained in references 12 and 13 as a result 
of this assumption. Balescu (ref. 14) uses a somewhat different method of 
linearizing the Fokker-Planck equation. He assumes that the test particles do 
not interact among themselves but only with the field particles, which are 
assumed to be in thermal equilibrium and, therefore, have a Maxwellian distri- 
bution. In Balescu's linearization, the distribution function of the test 
particles is not assumed to be close to equilibrium, however. Neither of the 
aforementioned linearization methods is used herein, inasmuch as the present 
analysis is intended to include those cases where the distribution function, 
although it is nearly spherically symmetric, is not close to a Maxwellian dis- 
tribution, and self-collisions are important. 

The simplified expansion method to be used herein proceeds in the follow- 
ing manner. The Lorentz approximation (eq. (27)) is substituted into the 
Fokker-Planck equation (eq. (21)), which is then reduced to two equations, 
describing the rates of change due to collisions of F u and F , respectively. 
In order to do this, it is necessary to determine at least the first few terms 
of the expansions for the functions H and G in spherical harmonics. It 
turns out that it is rather easy to expand functions H and G in full. 
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Rosenbluth, et al. (ref. 1), did this previously for the case of axial sym- 
metry. In the next section, the expansions of H and G in terms of general 
spherical harmonics are obtained by a method that is simpler than that of refer- 
ence 1. 


EXPANSION OF FUNCTIONS H AND G 

1 The function Hj.(v) defined in equation (9) involves the integration of 
g over the distribution function of the field particles according to 



where the angles 0 1 and cp T are shown in figure 1. In order to expand K^., 
it is convenient to first expand g" 1 as follows (ref. 8, p. 62, eq. (3.41): 


1 

g 


1 



z 



p^(cos r) 


(35) 


where v < (v > ) is the smaller (larger) of (v,v f ) and r is the angle between 

v and v', as shown in figure 1. The relation between and the angles 
9, 9, 0', 9' measured in the fixed coordinate system i, 1 k is 


cos T = v ! v ' 


w 


cos Q cos 0' + sin 0 sin 0 1 cos (9-9') 


(36) 


The addition theorem for spherical harmonics (ref. 8, p. 68, eq. (3.62)) is 
used to express P^(cos y) as a function of 9, 9, 9' , 9 1 


r) = 


£ 




*Zm 


lm K 


(37) 


When equations (23), (3b), and (37) are substituted into equation (34), it 
becomes 
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4jt 

21 * + 1 




2=0 ,m=- 2 


2 =0 m =-2 ' 



v 2 ' 7 

— ^ F^ m (v')Y z , ,(0,q>) 
2'+l f ^ 
v > 



sin 9' d0' Y* , , (0 ' ,cp ' )Y (0'/p') 
i m 2m 


(38) 


Because of the orthogonality property of the spherical harmonics (ref. 8, 
p. 65, eq.. (3.55)) 



d0 ’ " 5 7.7 & m. 


2'm 


2m' 


VI m’m 


(39) 


(where 5 , =1 if 2=2' and is zero if 2^2*), all terms in the m’ 

b b 

and 2' summations are seen to vanish when m' ^ m and V ) l, while the 
terms for m' = m and 2' = 2 integrate to unity. The resulting expansion 
for Ej. can he written 


00 2 


H t (v) 


■X v * “XXX 


s r (v) h ^ 


(40) 


f 2=0 m=-2 


where 




in m t l m f 

22 + 1 DLj. 


tf 




2+2 

— F 6m (v')dv' + 

Y-4+J- £ 


(v') 6 ' 1 ' ™2m 


/ 


— — ■ , -i F Zm (v' )dv’ 
(v 1 ) 2 " 1 f 


( 41 ) 
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Note that H^(v) is related to v 1 , 0, cp) by means of a one-dimensional linear 
integral operator in v 1 . 

The expansion of G^(v) is carried out hy first expanding g in Legendre 
polynomials. Both sides of equation (35) are multiplied "by 
g2 = v 2 + (v’)2 - 2w’ cos T, yielding 


g = 



v 


< 


Z+l 


v. 


[v 2 + (v') 2 ] 


P^ - 2w' cos Y 


> 


(42) 


If the recursion formula (ref. 8, p. 59, eq. (3.29)) 

(l + 1)P^ +1 - (22 + 1) cos r P z + ZP^i = 0 


(43) 


is used to eliminate cos Y P ^ from equation (42) and the numbering of the 
terms in the series is changed, the following expansion of g is immediately 
obtained^ 


g = 


00 _ 7 


i+i r 2 + <T '> 2 - 2TY ' (srrr ^ + rnh 


p^ ( cos r) 


2= 0 


(44) 


Since v^ + (v 1 )^ and w r are unchanged by an interchange of v and v', 
they can be replaced by (v< + v^) and v^v^, respectively, and equation (44) 
can be written more simply as 


g = 


00 

■ 

2=0 


v Z / V 2 \ 

< _ 22 - 1 < It, / 1 

l - ovt - 5 - P-,(cos r ) 




2 2 + 3 2 / 2 

>' 


(45) 


The argument used in deriving the expansion for Hj. is repeated for with 
equation (45) replacing equation (35). The result can be written down immedi- 
ately 


00 1 


G t ( ^ = Yj G f(v) = Y, Y, Yj G f m(v)Y Zm (0 ^) 

f “ 2=0 m=-2 


(46) 


This method of deriving the Legendre polynomial expansion of g was 
suggested by Willis H. Braun of Lewis. 
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where 



Note that G^(v), like H^(v), is related to F f (v T ,0,cp) h y means of a one- 
dimensional linear integral operator in v*. Equations (40), (41), (46), and 
(47) represent the extensions of equations (40), (41), (45), and (46) of refer- 
ence 1 to general spherical harmonics. 6 In the case m = 0, the two sets of 
results become identical. In reference 1, the expansions Hj. and are 

derived by first expanding g -1 and g in three-dimensional Fourier integrals. 
The subsequent integrations in that paper involve integrals of products of 
Bessel functions. The present derivation seems to be considerably simpler. In 
reference 3 the general terms of the expansions for Kj. and G^ are not given, 
so comparison is impossible. 

For subsequent purposes, only the terms for l - 0 and l = 1 in the ex- 
pansions of H{. and G^ (eqs. (40) and (46)) are required. These terms may be 

written in a form exactly analogous to equation (27) 

H t (v) = h£(v) + H^v) --+••• (48) 

G t (v) = G°(v) + Gj(v) * y + * • • (49) 


Since H^. and are related to F^(v',0,cp) by linear integral operators in 

v’ , it is clear that H®, G^?, and H^;, G^ are related to F^ and respec- 

tively, by the same integral operators, namely 




(v 1 "f" F®(v' )dv' 




)dv 


! 


( 50 ) 


®In equation (40) of reference 1, a^ a ^(v,n) should be replaced by 
A^ a ) (v,p) on the left side; a similar change should be made in equation (41). 
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^ (V) = ^^ (V) = ~ F tf ^2 / ( v ') 3 Pf(v')dv' + j 


->x 

F f (v , )dv' 


(51) 





\ 1 / — v ) — 
JT 


G?(v) = \ G°(v) = 4jt \ T 


tf 


v f (v') 2 F°(v')dv' + -L f 

Jo Jo 


V 


( v l )4j>0( yl ^ v l 


Z ! 


l 


+ / (v* )^ f9(v' ) dv' + 


if 


v 1 F°(v')dv’ 


3j(v) =") 3j(v) 


(52) 


Z 


4rt 

15 / Hf 


4r J (V ) 5 F^(v' )dv' 


5 / (v'^F^v'Jdv 1 

Jo 


+ v" 


/ CO ^.00 

F^(v ! )dY' -5 v / ( v ! ) 2 F^( v 1 )dv r 


(53) 


Equations (52) and (53) agree with equations (16) and (17) of reference 3, when 
equation (17) is written in vector form and the incorrect upper limit co in 
the first term of that equation is replaced h y v. No expressions for 
and K|. are given in reference 3. 


DERIVATION OF EQUATIONS FOR F° AND 


F 1 


Equation (21) will now he decomposed into equations for F and F „ The 
distribution functions F and F^ of test and field particles, and the poten- 
tials and G^_ are written 


F(v) = F°( v) + |(v) 
Fji( v) = F°(v) + |.(^) 


(22) 

(54) 


15 


I 



(55) 


Mv) = H$(v) + Tl(v) 


G t (v) = G°(v) + £(v) 


(56) 


where in the spirit of the Lorentz approximation S, i^ 9 i\ 9 £ are considered 

to he small perturbations. When products of £, i] 9 and £ are neglected 

in comparison with the first powers of these quantities 7 the various terms in 
equation (21) become 


F f F = F^F° + F^i + F u | f 


(57) 


aj L ^%_a^^f + a^3r 1 _ + ^ai_ 


(58) 


S 2 F 




■t 5 2 f° 5 2 g° a 2 F°_ a^ 


2^0 


b*G'. 




5v cTr dv Ov av dv. dv . dv . + SvTdv. Sv . av . + dv.dv. S77577 

i j ij ij i J ij t j ij t j 


(59) 


Since all the quantities with superscript 0 are functions of v, the follow- 
ing reductions can he made in the derivative terms: 


dF° _ dF° dy n dF° 

5v7 ~ 3v 5v^ ~ v Sv 


(60) 


a 

57 


2 f°_ _ a_ ('ll 8 f°\ 8 i.i SF° v i v .i a (i aF 0 ^ 

757 . 5v . \ v dv ) v 3v v S7 \v sw; 

i j J v 7 v 7 


(61) 


Similar expressions are of course obtained for H G and G G . If the gradient 
V y is written in spherical coordinates , it follows immediately that 


V i S77 


a -» „ a 

. = v * V v = v 37 


(62) 


Equation (62) is then used to show that 


FI 

1 3 


v i Y j 771777 = Y i 577 w 57" 


v. 5. . = v -7" 

i ij 377 dv 




V S .v 2 

37 ^2 

av 


(65) 


1 x j' J 

Expressions such as equations (60) and (61) are substituted for the derivatives 

0 
t 


of F°, H G , and C- G in equations (58) and (59), and use is made of equations 
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(62) and (63). The results maybe written 


dF dHf ^ ^0 ^ ^0 ^ ^ dg 

ov^ Sv^ civ civ dv dv dv civ 


( 64 ) 


d 2 F 


d 2 G + 


dv i dv j dv i dv j dv 2 dv 2 


d 2 F° S * G ? _2_ SG ? dF° 

v 2 37“ 


, 1 dF° 
+ yF 


^ £ + t a I + V a(i^?)a!i (6S) 

V b dv \v dv J ^ y 2 v dv v dv \v dv J ^2 v ' 


If the quantities |, t), £ are identified as the terms for 2 = 1 in the 
spherical harmonics expansion of F, Hp, and Ch , it follows from equations 
(27), (48), and (49) that 


!(v) = F G (v) • £ 


ri(v) = H f (v) • ^ 

C(v) = fli(v) * | 


( 66 ) 

(67) 

( 68 ) 


The derivatives of £ , t^, £ with respect to v can be calculated immediately, 
since the unit vector v/v remains fixed in magnitude and direction as v 
changes but 0 and cp remain fixed. Thus, for example 

as = v 

Sv V 


dv* 


dF^v) 

dv 

(69) 

d 2 F G ( v) 

(70) 

dv 2 

7 2 | and V 2 
v v 

£ is to note that, for 


* = 



^Zm( v )Y ( 0, cp) 

Im 


(71) 


2=0 m =-Z 


the Laplacian is 
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7=0 m=-Z 


2 ^ 
civ J 


r 1 a / . a \ 1 a 2 1 

2 sin e ¥ ( sin 9 'be) . 2 . 2 I 

v [_ ' / sin 0 ocp J 


’7m| 


But Y satisfies the equation (ref. 15, p. 53, eq. (5)) 


[ sl " 6 a ® ( sln 15 *>) + sin 2 e a, 2 ] 


+ -Vqjl r im = - !(l + 1)J ln 


hence y 


V 2 1)/ = 

Y Y 


o° p _ 7- 

Z=0 m=-Z 


^Zmy 


Zm 


Then, "because | can "be written (see eq. (26)) 

1 

| = £ F lm (v)Y lm (e, 9 ) 


m=-l 


it follows from equation (74) that 


V 2 i = 


F lm( v)Yi m (e,cp) 


If (- 2 If) - 7] 1 

L J m=-l 

[M (■>*)- -*][«*>-;] 

b s (•* i) - >]»><” 


= I.li 

V 


(72) 


(73) 


(74) 


(75) 


(76) 


The expression for V 2 £ is similar to equation (76), with Fl replaced by 

G^. Substitution of equations (69), (70), and (76) and their counterparts for 
the derivatives of tj and £ into equations (64) and (65) yields after some 
rearrangement and cancellation of terms 
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3? SH f 


aH f + 1 

dv dv v 


( 


3f° 

5v~ dv 


+ aa ? ^ 

dv dv ) 


(77) 


S 2 f a2s t 

dvj_3vj dv-^dvj 


a 2 ? 0 a 2 a° ^ 2 3 °° a ? 0 ( v 

£ v 2 £ v 2 v 2 ^v" v 


_2_ dF° /^ G t 1 

v 2 dv \ dv~ v t / 


_2_ ^ G t ( 3F 1 _ 1 pl\ + d 2 F° + ^ 2Q t d 2 ! 1 

v 2 \ dv v / dv^ dv^ dv^ dv^ 


(78) 


When equations (57), (77), and (78) are substituted into equation (2l), the 
Fokker-Planck equation can finally he decomposed into the following two compo- 
nent equations for F^ and F^: 


(i 2 ) -3 

x /cc -p L 


™t n ^O-oO ^ dF° / m f " ™t ^ H f . l ^ G f 


4n -1 F tf FfF^ + [ — - ^ + — _ 

mf UJ - x dv l m-f + m^ dv v 2 dv 


+ 1 ^^^f d 2 F^ 
^ dv 2 dv 2 


- 

) 


(79) 


/a? 1 ) 

C' i ^ G f d 2 ]? 1 

dF 1 1 

f m f - m t dH G x oG G \ 

V^cc ‘ , 

/ > _ 2 dv 2 dv 2 

+ 3v~' 

\m.f + m-f dv v 2 dv / 


+ F- 


/V ^ 0 l SG f ^ I l d 2 F° 

'V m f P tf F f ' v 3 ) 2 Sv 2 Sv 2 


. dF° (] 


'mf 


m . 


(80) 


COMPARISON WITH OTHER WORK 

Equations (79) and (80) are not exactly in the form given by Dreicer 
(compare eqs. (18) and (19) of ref. 3). Dreicer evidently eliminated the poten- 
tial H by means of relations (17) and (18). To compare the results obtained 
herein with Dreicer T s, it is necessary to do the same thing. From equation (17) 


H* - 


m^- + mf 
ZXQf 


o f 


(81) 
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The expansions for (eq. (40)) and for Gf (eq. (46)) are substituted into 
equation (81), and the Laplacian of Gf is calculated. When equation (74) is 
used and the coefficients of Y^ m on both sides of equation (8l) are equated, 
the result is 



Equation (83) is substituted for in equation (79), -which, after some rear- 

rangement, can be written as 




4«br tf J?0]?o + §£ 


\v2 


dG° 


+ ~ 


2c d G f 


Sv 4 


+ c 


5 3 g° 


) + l 


c^Gf ^2jp0 


(85) 


where 


b = 



( 86 ) 


c 


2mf 


(87) 


Equation (85) is exactly the same as Dreicer’s equation (18). It is not clear 
why he chooses to eliminate Hf in favor of Gf, since higher derivatives are 
obtained in this manner. 


In order to compare equation (80 ) with Dreicer’s equation (19), the follow- 
ing relation must be used: 


- 8jtr tf F f (v) = vj Gf (88) 

which follows from equation (18). When the expansions of F f and G f (eqs. 

(23) and (46), respectively) are substituted in equation (88), the Laplacian 
operation performed on Gf twice, with the use of equation (74) each time, and 
finally the coefficients of Yj m on each side of equation (88) are set equal to 
each other; the following result is obtained: 
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- 8«r tf F^(v) . £ (v 2 of (v) (89) 

In particular, the relations that are needed are 

- 8*r tf F!9(v) = [i |L(v 2 ^ G°(v) (90) 

- 8«r tf ^(v) = [-§ (v Z Jj) - 3 2 (v) (91) 

Equations (90) and (91) are substituted for F^ and F^; in equation (80). 
Also, and H^; are eliminated by means of equations (85) and (84). After 

some manipulation, equation (80) becomes 


(&\ = V h 

VSt / / - 2 


S G^ £)2]?1 

3v 2 Sv 2 


Sf 1 / 

Sv v 


s 3 g« 2c a 2 o° 

C r — + — + 

c)v 3 V V 


± d 3.) 

V 2 / 


->-i /b 

S 4 G? 

2b 

S 3 G? 

1 

SG?\ 

i S 2 f° 

a 2 ^ 


- ^2 

Sv 4 

V 

Sv 3 

V 3 

377 

' 2 Sv 2 

Sv 2 


Sf° i 

( 8^ 


2c 

2b - 1 

dGf i 

- 2b 

7?l\ 

+ 3T 1 

^Sv 3 

+ ■ 

v Sv 2 

+ - 

V 2 

• 37” + - 

V 3 


- 

S 4 gJ 
^ + 

2b 


2b 

S 2 C* N 

)] 

(92) 

\2 

Sv 4 

V 

Sv 3 

V 2 

dv 2 y 

'J 



which is in the form of equation (19) of reference 3. It would appear that 
equation (80) is more convenient than equation (92), since the latter contains 
higher derivatives of the G functions. 

Comparison of equation (92) with Dreicer’s equation (19) discloses two 
differences. First, and more important, the last term in Dreicer's equation, 

-(4/v 3 )dG^/£v, was not obtained herein. Since the calculations were checked 

by a number of different methods and the term in question was unobtainable, it 
is believed to be incorrect. The second difference is evidently a misprint in 

the term (l/2) (S 2 F°/Sv 2 ) (S 2 (J^/Sv 2 ) , where Dreicer has F^ instead of F®. 
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Shkarofsky (ref. 4=) quotes Fokker-Planck collision terms similar to equa- 
tions (79) and (80), but specialized to the case of electron- electron colli- 
sions. A comparison of the present results with those of reference 4 in that 

special case showed exact agreement for (33^/dt) (ref. 4, eq. (5b)), but 

wide differences in the expression for (5F°/dt) ee (ref. 4, eq. (5a)). 

Shkarofsky’ s equation (5a) is in fact dimensionally incorrect, and the term con- 
taining the second derivative of F^ is missing from it. 

Finally, some comments on an early paper by Cahn (ref. 5) may be of inter- 
est, because he apparently managed to obtain an electron- electron collision 
term in a form that enabled him to integrate the Boltzmann equation and express 
the distribution function by means of a single integral. If this result were 
correct, it would be quite significant due to its simplicity; however, it is 
believed that Cahn l s result is not correct. The reason has to do with the way 
in which Cahn carried out integrations over the velocities v of the field 
particles. In reference 5 (p. 298), the lower limit of integration on the 
relative speed [v* - v| is quoted as v* - v. Actually this should be 
jv 1 - v|; otherwise, it is incorrect in the case where v* < v, because the 
relative speed must be positive. The result of the angular integrations is 
therefore different, depending on whether v 1 < v or v* > v. The integration 
over v 1 from 0 to 00 then breaks into two integrals with different inte- 
grands, one integral from 0 to v and the other from v to °o, as obtained 
in this report (see eqs. (4l) and (47)) and references 1 and 3. Cahn has only 
one integral, from 0 to For instance, the following result is given in 

reference 5 (p. 299): 



where Cahn T s notation need not be explained in detail here, except to note that 
£ is a constant related to F , while bPO involves the distribution func- 
tions of the field and test particles and their derivatives. The result should 
be 



This difference is important, because Cahn l s electron- electron collision term 
is the divergence of the vector J±, and a number of extra terms would appear 
from the differentiation of the integrals having variable limits. These extra 
terms complicate the result considerably, and there is no longer a possibility 
of obtaining the distribution function F° by simple quadratures, as claimed 
by Cahn. 


22 


iihiii 1 11 11 11 ilia 1 ail 111 


CONCLUDING REMARKS 


The first two terms in the spherical harmonics expansion of the Fokker- 
Planck equation have "been obtained. The results are nonlinear integro- 
differential equations whose variables are the speed v of the test particle 
and the time t. Although neither these results nor the general method of de- 
riving them is entirely new ; the fairly complete and straightforward deriva- 
tions presented herein resolves the errors and differences that exist in previ- 
ously published results. Although the Fokker-Planck equation is the best pres- 
ently available means of including the effect of electron-electron interactions 
on the distribution function, it is still rather intractable to all but numeri- 
cal means. A simpler electron-electron collision model, which would bear the 
same relation to the Fokker-Planck collision term as the Krook model does to 
the Boltzmann collision integral, would be invaluable for workers in the field 
of plasma physics. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, January 8, 1965. 


23 



APPENDIX A 


SYMBOLS 

The niks system of units is used throughout this report 
i'k* 1 component of acceleration 
B magnetic field 

b m t /m f 

b Q impact parameter for 90° scattering, eq. (14) 

c (m-f - mp /2nif 

d Debye radius, eq. (13) 

E electric field 

e electron charge (e = 1.602x10"-^ C) 

ej» electric charge of field particle, ±Z^e 

e^ electric charge of test particle, ±&j.e 

F distribution function of test particles 

F|. distribution function of field particles 

G Hosenbluth potential, defined in eq. (10) 

g relative speed of test and field particles, eq. (3) 

— > — > -> 

g relative velocity, v - v 1 

g 1J conjugate of metric tensor, eq_. (B5) 

gj_j metric tensor 

H Roseribluth potential, defined in eq.. (9) 

Christoff el symbol of second kind, eq.. (B2) 

Im[ ] imaginary part of 

AAA 

i,j,k unit vectors along Cartesian coordinate axes 
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[ij,k] 

k 


DLf 


m. 


m 


tf 


N 




Re[ ] 


r 

s 

T 

t 

u 1 

V 

V* 

-> 

V 

V* 


Av ± 


v! 


x i 


-Zm 


Cliristoffel symbol of first kind, eq. (B3) 

Boltzmann constant, 1. 3804xl0~^ 3 j/°K 
mass of field particle 
mass of test particle 
reduced mass, m^_m^/ ( m^_ + m^) 
particle number density, m"3 
Legendre polynomial of order l 

associated Legendre polynomial, defined by eq. (24b) 

real part of 

position vector 

arc length in velocity space 

temperature 

time 

general curvilinear coordinate in velocity space 

magnitude of velocity of test particle, (v^v^ )l/2 

magnitude of velocity of field particle, ( v_j, v_p A 

velocity vector of test particle 

velocity of field particle 

drift velocity of test particles 

i*kk component of velocity of test particle 

increment in v^ due to collision 

i^* 1 component of velocity of field particle 

larger of (v,v ! ) 

smaller of (v,v T ) 

i^ n component of position vector r 

spherical harmonic function defined by eq. (24a) 
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x tf 

r 

8 . ., 8 * 

J 


5 

T1 

e 

^tf 

M- 

l 

i 


f 

O' 

9 

X 

a 

( ) 


V 

V 2 

v 


charge number 

defined by eq_. (ll) 

— ^ ^ 

angle between v and v* 

Kronecker deltas in Cartesian and general coordinates , respectively 

dielectric permittivity of vacuum, 8.854xl0“ 12 P/m 

perturbation in eq. (56) 

perturbation in eq. (55) 

polar angle in velocity space 

defined by eq. (12) 

cos 6 

perturbation in F, eq, (22) 
perturbation in F , eq. (54) 


differential scattering cross section 
azimuthal angle in velocity space 
scattering angle in center of mass coordinates 
scattering solid angle 
average value 

gradient operator in velocity space 
Laplacian operator in velocity space 
Subscripts : 
c collision 

cc Coulomb collision 

ee electron-electron collision 

f field particle 

h,i,j,k tensor indices of covariant vectors and tensors 

s species s 
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t test particle 

, denotes covariant differentiation 

Superscripts : 

h,i,j,k tensor indices of contravariant vectors and tensors 

2 index denoting term in Legendre polynomial expansion 

2m indices denoting term in spherical harmonics expansion 

m index denoting term in azimuthal angle expansion 

0 zero-order terms 

1 first-order terms 

** complex conjugate 

-> vector 

unit vector 
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APPENDIX B 


TRANSFORMATION OF THE FOKKER-PLANCK EQUATION TO 
SPHERICAL COORDINATES 

In this appendix, equation (21) is transformed to the spherical coordi- 
nates (v, 0, cp) in velocity space. This transformation was performed on the more 
complicated equation (19) in reference 1 for the special case of axial symmetry. 
Basically the same method is used, namely, transformation to a general curvi- 
linear coordinate system and the specialization of the result to spherical co- 
ordinates. Axial symmetry is not assumed, however. 

Let the general curvilinear coordinates he denoted as u- 1 -, u 2 , u 3 . First 
the term (d 2 F/dv.[dvj ) , which was written as a Cartesian tensor in equation (21), 

will he transformed into a second-order covariant tensor. Since F is a 
scalar, the proper extension of (d 2 F/<dv\j_dVj) to general coordinates is the sec- 
ond covariant derivative of F (ref. 16, p. 32, eq. (22.7)); that is, 


' ^ Su' 


fF f hi 5F_ 

' l 1 ')/ 3u h 


•where (ref. 16, p. 26) 


&}- 


g^[ ij,k] 


and 


r .. vl Sg i 3 


(Bl) 


(B2) 


(B3) 


are the Christoffel symbols of the second and first kinds, respectively. The 
quantities g * a are the elements of the metric tensor; that is, the element of 

1 J “I p 'Z 

distance ds between two points whose coordinates differ by du x , du^, du° is 

ds 2 = g. .du^du^ (B4) 

and g^ is the contravariant tensor conjugate to g. ., defined such that 

g«g. k = 8j (B5) 


where 6^ is the Kronecker delta. 
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In Cartesian coordinates. 


(S 2 G t /Sv ± Svj) 


is defined such that 


<Av . Ay. ) 


d v . dv . 
i 0 


(B6) 


hut, in general curvilinear coordinates, this^must he written as a second-order 
contravariant tensor, because the velocity u^ and velocity increment AuA 
are contravariant vectors. Thus the proper extension of equation (B6) to curvi- 
linear coordinates is 


(Au^-AuJ) = gihgjkQ. (B7) 

where the nontensor subscript t on G should not he confused with the tensor 
indices h,i,j,k. 

The term (Sf/Sv^) (bE^/bv/) in equation (21) can he transformed easily, 
since the first covariant derivative of F is simply 


F . 

J 1 


5f 

Su 1 


(B8) 


and the quantity (Av^) = bE^/bv^ defined in equation (7), when written as a 
contravariant tensor, becomes 

(Au 1 ) = g^H (B9) 


In general curvilinear coordinates, equation (21) becomes 


(i) oc - 


^ r tf P f (v)F(?) 


- ^ 

8 ctt : 


'AZ 


m t - m f i 


m 


t % + 2 s ih s JkF ,ij G t,hk ( B1 °) 


For spherical coordinates the element of distance ds is given by 


ds 2 = dv 2 + v 2 d 0 2 + v 2 sin 2 0 dcp 2 = dv 2 + — X djj 2 + v 2 (l - p. 2 )dcp 2 

1 - ia 2 


(Bll) 


where for convenience the coordinate 
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|JL = cos G 


(B12) 


will "be used in place of 0, as in references 1 and 3. Hence , if 

u-'- = v, u^ = p, = cp 

it follows from equations (B4), (B5 ) , and (Bll) that 


§ 11 " 11 " 1 


S 22 “ 22 


v 2 


1 - p 


S 33 33 v2 ( 1 " ^ 

§ 


L , = g lJ = 0 (i f J) 


(B13) 

(B14a) 

(B14b) 

(B14c) 

(B14d) 


Equation (B14d), which is true for any orthogonal coordinates, mayhe used to 
considerably simplify expressions (B2) and (B3) for the Christoffel symbols 
(see ref. 16, p. 28, Ex. 2). For spherical coordinates, it then follows from 
equations (B2), (B3), and (B14) that 


fe) ' ' 7T7 <B15a) 

{22} ■ - {25} - - {I2} = 77^2 < B15b > 

= - ^(1 - p 2 ) (B15c) 

| 3 2 3 | = p(l - p 2 ) (B15d) 

{-} - h - {-} = I-} - * (— ) 


all other Christoffel symbols being zero. Substitution of equations (B15) into 
equation (Bl) yields 


F 



(B16a) 
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(B16b) 


j> _ 5 F , v SF (j. dF 

,22 ~ . 2 . 2 5?" 2 51 1 

dpi 1 - pi 1 - |-i 


,33 - 2 

' dcp 


5 2 f 2 v 5f o. 5 f 

+ v(l - M-0 ^ - M-U - M-0 


5T 


F ,12 F ,21 


5 2 F _ l 5 f 
dvdq v 55 


F = F 
,13 ,31 


5 2 f 1 5f 

Svckp v Sep 


2 

_ 5 F q SF 

,23 ” ^,32 ” 5^9 |^2 ckp 


F _ = F 


The expressions for the components of G. . . are of course analogous to 

D , 1 J 

formulas of equation (B16). Substituting equations (B14) and (B16) into 
tion (BIO) results in 


(»L ■ X- 


m+ m.p - mi 

4tc — I\ p (v)F(v) + — 

m_p tf f w v ; m , + m„ 
f t f 


dr 9% 1 -Jig dF SH f 1 dF SH f 

dv - _2 5jT <5~ + dcp dcp~ 


1 dfT d G f + (1 - u 2 ) 2 / ?pF + v dF _ n 3 f\ ( 5 G f + v 5G f y_ 


2 dv 2 dv 2 ' 2v 4 W 2 ' 1 - ^ ^ ' 1 - n 2 ^)W 1 - M 2 ^ 


1 - l 


2v 4 (1 - pi 2 ) 


d 2 F 


dcp 1 


.ft 2 ^ dF , , 2\ dF 

j, + tU - n ) ^ - n(l - li ) gj;- 


d 2 G.p 


dcp d 


O SG f 

+ v(l - ^ 2 ) - pi(l 


:Lzjd(^F 1 dFV ^ 2 °f 1 d&A 1 / d 2 F 1 1 

v 2 \dvdpi v dpi J ydvdpi v dpi J + _ ^ 2 ) ydvdcp " v depi ydvdcp v dcp 


, 1 fd 2 F u drV^f n 5G f 

v 4 \dpidcp 1 _ ^2 dcpyydpid 9 1 _ ^2 dcp 


(B16c) 


(B16d) 

(B16e) 


(B16f) 

the 

equa- 


_ SG d 
2 ^ / 



(B17) 
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Equation (B17) is the Fokker-Planck equation in spherical coordinates. It 
is possible to utilize it to obtain the general term in the spherical harmonics 
expansion of (SF/St) cc . This expansion could in principle be obtained by sub- 
stituting the following expansion in equation (B17): 

F(v) = £ V F Zm (v)Y (0,cp) (23) 

Z=0 m=iz 

together with the analogous expansions for H^, and G^. The recurrence 

relations for the functions Y^ (ref. 17, p. 181) would he utilized to express 
their derivatives and the products of their derivatives with the quantities 
p, (1 - p 2 ) 1 / 2 , and so forth, in terms of the functions Yj ffl themselves. It 
would then he necessary to expand products of spherical harmonics as follows: 


Y ( 0, <p) Y (0,cp) 

Vl J 2 n 2 


00 



z=o 


(B18) 


where the constants 


C? m 7 m 7 are related to the Clehsch-Gordon coefficients 
T12F 


(ref, 17, p. 514). Finally, the result would have to be rearranged so as to 
separate out the various terms in Y^ m . Although the procedure is straight- 
forward, the amount of calculation would be rather large, and the usefulness of 
the general term in the expansion would probably not justify the labor in ob- 
taining it. For this reason the present report is restricted to the first two 
terms (Z = 0,1) in the spherical harmonics expansion of the Fokker-Planck equa- 
tion. The nonlinearity of equation (B17) is responsible for the difficulty in 
working out the general term of the expansion. The expansion of a much simpler 
linear collision term appropriate for collisions with infinitely massive field 
particles (the so-called Lorentz model) has been carried out in reference 18. 
Reference 18 also contains the spherical harmonics expansion of the electric 
and magnetic field terms and the spatial gradient terms in the Boltzmann equa- 
tion. 
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